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Abstract 



We calculate the decay rate for inclusive B tvX decays in a two Higgs 
doublet model using heavy quark expansion and operator product expansion. 
Combined with the recent measurement of Br(i? tvX), we find a limit 
tan/3 < 0.58m///GeV at 90% C.L.. 
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One of the most exciting topic in the Standard Model (SM) is the existence of the Higgs 
boson and its nature. Even though there is only one neutral Higgs boson in the SM, more 
complicated structure of the Higgs boson sector is usually assumed in many models beyond 
the SM [0. The two Higgs doublet model is one of the minimal extensions of the SM. One 
or more charged Higgs bosons inevitably exist when the multi-Higgs doublets are assumed. 
Hence it is interesting to find any of the charged Higgs boson effects in probing new physics. 
With the CLEO bound, the flavour- changing process 6 — 57 is known to provide the most 
stringent constraint on the charged Higgs parameter space 0. It is also reported recently 
that the vertex Z —>■ bb with one - loop corrections can also provide additional constraints 
for tan/5 < 1 0. 

The inclusive semileptonic decays of B meson involving a tau lepton in the final state can 
be also an useful probe for the charged Higgs effects of multi-Higgs models. This process is 
one of the latest measured channel which constitutes a significant fraction of all semileptonic 
B decays. The ALEPH and L3 collaborations at LEP e^e~ collider reported on the improved 
measurement of branching ratios of the decay B tvX |^,^ : 

Bi{B Xtu) = (2.76 ± 0A7{stat.) ± 0A3{syst.))% ALEPH 
Br{B Xtu) = (2.4 ± 0.7{stat.) ± 0.8{syst.))% L3 

Being distinct from other probes of charged Higgs effects, this process is not sensitive to 
the top quark mass. Furthermore we can do the systematic study for the nonperturbative 
contributions of this decay process in the framework of heavy quark effective theory (HQET). 

Recently, there has been a great progress in systematic study of inclusive semileptonic 
decays of hadrons containing a single heavy quark. In the pioneering work by Chay, Georgi 
and Grinstein 0, the method by means of the operator product expansion (OPE) in the 
heavy quark effective theory was suggested, which is motivated by deep inelastic scattering. 
Following this formalism, one can implement the QCD - based expansion in powers of 1/ mg 
to calculate inclusive B WX decay rates and lepton spectra. The leading order terms 
coincides with those of the free quark decay. The noticeable result is that the leading 

2 



nonperturbative corrections are of order of l/mq because terms of order of l/mq vanish. 
Several authors have studied the inclusive B —>■ tvX decay rates in the framework of the 



SM 10.11.12 



Here we calculate the decay rates and lepton spectra for the the inclusive semileptonic 
decays of B meson into a tau lepton with the help of the HQET in a two Higgs doublet 
model, so-called type II. In the type II model, one Higgs doublet couples to up type quark 
sector and the other Higgs doublet to down type quark sector, while only one Higgs doublet 
couples to the fermions in type I model. We analyse only type II model because it is required 
in the minimal supersymmetric standard model (MSSM) and type I is less interesting as 
type I is similar to the SM. We also discuss their phenomenological meaning as constraint 
for the charged Higgs boson. 

In the two Higgs doublet model of type II, the inclusive B tuX decay is due to the 
hamiltonian density, 

Ti = Tiw + T^H (1) 



where 



'Hw — — ^^^^^ ^M'^T 



Hh = -^^V,t JJr (2) 



H 



with 



= c7m(1 - 75)^ 
J = c KX(1 + 75) + mcY{l - 75)] b 

Jr = f [mrX{l - 75)] Ur (3) 



The ratio of two vacuum expectation values are defined by 



X = tan/5 = — 

Vl 



and Y — 1/X. The inclusive differential decay rate for the process B{pb) — > 

xdq^ dEr dE^ (4) 
is related to the three hadronic quantaties, 

W^'' = Y.{2t:)H\pb - q-px){B{v)\J^'^\X){X\r\B{v)) 

X 

W = Y.{2ti)H\pb - q-px){B{v)\r^\X){X\J\B{v)) + H.c. 

X 

W = Y.i'^nf^'iPB - q-px){B{v)\j'^\X){X\J\B{v)) (5) 

X 

where leptonic tensor L^^ is given by 

L'^" = PrV + P/Pu'' - g^'Pr ■ P9 + le^^'^yrpPu, (6) 

and q = Pr + Pp \s the momentum transfer. We can expand the quantities in terms of form 
factors. The hadronic tensor W^^ is parameterized by 

W"" = -g'^'Wi + v''v''W2 - ie^^'^'^VaqpWs + q^q^W^ + {q^'v'' + q''v^')W5 , (7) 

and 

W^" = qf'Vi + v^V2 , (8) 

and 

W^S. (9) 

where v is the four velocity of B meson. 

Now our goal is to calculate the quantities W^", W^, W . The form factors in W^", W^, 
W are determined by following amplitudes T'^'^, T'^, T respectively, 

Tf"" = -i j d^xe-"i-''{B{v)\T{J>'\x)r{Q)}\B{v)) 
T^'^-iJ d^xe-''^-^B{v)\T{J^'\x)J{0) + H.c.}\B{v)) 
T=-i j d^xe-"i-''{B{v)\T{j\x)J{Q)}\B{v)). (10) 
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The time ordered products of two currents 

-i J d^xe-'^-''T{JJ] (11) 

can be computed using an OPE in terms of operators involving b-quark fields in the HQET. 
At tree level, one may obtain OPE up to operators of dimension five by considering the rel- 
evant diagrams in fig.l. The calculation is straightforward and discribed in detail elsewhere 
Pfljlji]- Here we just present the results. Introducing a set of dimensionless variables 



.2 2Ey 



s2 _ 

rrii,^ "nib ''^6 



and defining the mass ratios 

Pc = — ^, pT = — J , (13) 

where Ej. and are the energies of tau lepton and neutrino in B rest frame, the triple 
differential decay rate can be written by 



Vbdq^dydx \ y + ^Jy'^ - Ap^ ) \y + ^Jy'^ - Ap^ J 

{ {e-Pr)m + ^{yx-f + Pr)W, 

+ l[fiy -X)-Priy + X)]W3 + ^Prif " Pr)Wi + PrXW^ 

-lR^/p^[iq^-pr)Vl + xV2] 

+\R\f - Pr)S } (14) 



where 



r, 



b 



1927r3 
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R = . (15) 

The expressions for the form factors Wi, Vi, and 5* at tree level and to order 1/mf in the 
operator product expansion are given in the appendix A. We assumed X > 1 and neglected 
the terms propertional to Y^p^. The leading term in l/rrib expansion reproduces the parton 



model results. The nonperturbative corrections of order is written in terms of two 

hadronic parameters Ai, A2. The parameters are related to the kinetic energy Kh of the 
b-quark inside the B-meson and to the mass splitting between B and B* mesons. 

Integrating over the kinematic variable and x, we obtain the differential decay rates 
with respect to the rescaled lepton energy y which can be written as 

= Ji-RJ2 + R^J^ (16) 

where the rather lengthy expressions of Ji^ J2 and are given in the appendix B. 
Finally, we perform the y-integration over the kinametic region 

< y < 1 + Pr - Pc (17) 
to obtain the total rates. For F, we find 

^ = Fi - i? F2 + i?'F3 (18) 
J- h 

where again the expression of Fi, F2 and F3 are given in the appendix C. As can be seen 



easily, our result reduces to the SM result of the ref. [10,11,T2[ if i?=0 



From the above result, the following relations holds 

Y{B^rvX) = Y'^^^\B ^TvX){\- AR^BE^) (19) 

where A = F2/F1 and B = F3/F1. Using this relation, one can obtain the constraint 
on tan (3 /rriH from the measurement of branching ratio, Br (5 rvX ). To obtain the 
numerical value, we need the values of the tau mass m^-, the heavy quark masses m;, and 
mc, the hadronic parameters Ai and A2 and CKM matrix element \Vcb\ as input parameters. 
From the heavy quark expansion of B and D meson masses, four parameters m?,, mc, Ai, A2 
are reduced to three parameters A, Ai, A2, where A is associated with the effective mass of 
light degrees of freedom in heavy mesons. We then take the values of them used in the paper 
of Falk et al. [110. To avoid the uncertainty in CKM matrix element \Vcb\i we normalize the 



branching ratio oi B tvX to that of S ^ evX ||T3| : 
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Br(E ^ ei/X) = (10.7 ±0.5)%. (20) 

including the perturbative QCD correction. 

This procedure leads to the following values of A and B: 

A = 0.43, B = 0.27 (21) 

With the measured branching ratio of L3 |^ , this gives the refined constraint 

A- < ^ (22) 

at 90% C.L.. 



Isidori [|T4| has estimated this limit with free quark decay formulation and the first 



measurement of this channel in ALEPH [15|. The value obtained by them was X < 



O.bAmH /GeV. But they made a mistake in numerical calculation of A in eq. (25). With 
correcting the value of A, their result should be modified as X < 0.61m///GeV. The dif- 
ference of values between their result and our one is principally caused by the fact that we 
used the recent data to estimate our limit. And it is also due to the differences of values of 
input parameters as well as l/m^, corrections. 

In conclusion, we calculate the decay rate of inclusive B — * tuX decayin two Higgs 
doublet model using HQET and OPE. From the recent measurement of Br( B tvX), we 
get the limit on tan jS/rriH. 
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Appendix A 



The form factors Wj, t^, and S to order 1/ml in the operator product expansion are 
given by 
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#. = 5(I){1(2-,-.)-^, 



[8f + 6(y + x)-5(y + xf] 



+ + 14(y + x) - 5(y + xf - 16]} 



'■6 



5(Ai + 3A 



6mi 



(. + .) + ^[5(. + .)-4]} 



2 



+^'(^^^12;;? (. + .) + 4[5(y + .)-12]} 



2(Ai + 3A2) 



^1 r./ X 01 5A: 



+8] + ^ (.+.)} 



A 3A 

+^'^^^^6;^ [-^ " + + 24 [2 - + 
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S = S{z){[l -l{y + x)] + [4 - (y + x)]} 

+(?/+-) -^(y+^ni 




where 

i = l + g^-Pc-2/-a; + ie. (A. 2) 
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Appendix B 



The expressions of the lepton energy spectra J2 and Jz are given by 



X 



[2/^ - 3y(l + p.) + 8p.] + xi[-?,y' + 61/(1 + pr) - I2p,] 



[bxi{y' - V(l + Pr) + 2y(3 + 7p,) + 4p,(p, - 5)) 



+ pr-y) 

+3xo(-5|/^ + 1/2(17 + 15p,) - 2/(24 + 46p,) - p,(18p, - 70)) 
+3(5y^ - 102/2(1 + Pr) + 42/(3 + 4p0 + 16p,(p, - 2))] 



Ai 



,[3(2/^ - 2y\l + p,) + %ypr{l + p.) - IGp.^^ 



3m^(l + p^ -2/)2' 
+2xg(y^ - 5y=^(l + pr) + 2y2(5 + lip, + Sp,^) - 40yp,(l + p,) 

-2p,(5-38p, + 5p,2)) 

+3x2(-22/^ + 82/^(1 + p,) - 2/2(15 + 28p, + ISp,^) + 52yp,(l + p,) 
+18p,(l-6p, + p,2)) 

+6xo(y^ - 3y^(l + pr) + y\b + 6p, + Sp,^) _ i2yp^(l + p^) 

-8p,(l-4p, + p,2))]j , 



J2 



:^y^-Apr 



xl (12 + 12p,-122/) 
A2 



[a;o(-362/2 + 362/(3 + pr) - 72(1 + p,)) 



ml{l + pr-y) 

+a;2(92/2 - 362/ + 27 + 18p, - 9p,2)] 

Ai 



ml{l+ pr-y) 



-XAy' - 42/2(1 + Pr) - 16ypr + 16p,(l + p,) 



+XoiAy^ - 162/2(1 + Pr) + 22/(6 + 28p, + 6p,2) _ 32p^(i + p^)) 

+2^o(y' - + Pr) + y(3 - lOp, + 3p,2) 

-3 + 7p, + 7p,2 _ 3^^3)|1 ^ 



xl (2/2 - 32/(1 + Pr) + 8p,) - 3a;2(2/2 - 3y{l + p,) + 4p,) 
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with 



A, 



2ml{l + pr-y) 



+3x2(13|/3 - y\Al + 39p^) + 2y(15 + 64p^ + 3^^^^ - 4p^(19 + 2>pr)) 
-xl{W - 521/2(1 + pr) + 1/(69 + 182p. + 9p.2) + i6p,(14 - p,)) 

'^l r r, / .4 o„,3/i I „ ^ I o„.„ (-1 I „ ^ ii5„ 2^ 



?,ml{l + Pr - y) 



3 (y^ - 2t/3(l + p,) + %ypr{l + p,) - 16p,2) 



+3xo(2?/^ - 6i/3(l + p^) + 2|/2(5 + 9p, + hpr^) - Qyprib + 7pr) - 4p,(4 - 19p, + p,^)) 
+3xl{-Uy'' + y\S3 + 56p^) - y'(123 + 304p^ + 75pr^) 
+1/(45 + 358p^ + 361p^2 + 24p^3) _ Qp^(^ii + 50p^ + ISp^^)) 

+2xl{7y^ - 35y%l + pr) + 2y^{2Q + 77p^ + 26p^2) 



-22/(9 + 95p, + 95p,2 + 9p,=^) + 2p,(19 + 86p, + Wp,^))] 



(B.l) 



1 + Pr - y 



(B.2) 
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Appendix C 



The decay rates Fi, r2 and Fa are given by 



Fi = \/A 



(1 + - 7(Pc + Pr) - + Pr^) + (P/ + Pr^) 

+p,p^[l2-7(pe + Pr)]} 

3A 

f;^{-3 + 5(pc + Pr) - l%Pc + Pr') + 5(p/ + p^') 

+ 7pePr[4-5(pe + Pr)]} 



+12(1 + 

-12(1 + 
F2 = v^a/A 



Ai + 3A2 
2ml 



)[Pc'ln 



(1 + Pe - Pr + VA)2 , 2^ (1 + p^_p^ + VA) 



4pc 



+ Pr In- 



4pr 



-)pcVr'ln 



(1 - Pr - Pc + VA)^ 
2 - lOpe + 20pr - 4p,2 + 2pr' - lOp^Pr 



Ai + 15A 2^ ^ 2 „ 2i 
2m? 



■f-^(-l + 5pc - lOpr + 2pc^ - Pr' + ^PcPt) 

3A 

+ :r4(-9 + 21Pc - 30pr - 6pe' + 3pr' " 15pcPr) 



Ai 



3A2 



+3{pe'[4(l - p,) - ^(1 - p,) - ^(1 + 3pr)}ln 



(1 + Pe - Pr + 



mi 



^Pc 



2Ai 



+3{-8pr + 16pcPr - 8pc Pr " 8pr H jip-r ~ ^PcPr + Pc Pr + Pr ) 



ml 

, 6A2._ „ Q 2 , 2.1 (1 + Pr - Pc + a/A)^ 
] 2"(5pr - 2pcPr - 3pc Pr + Pr jm 

ml 4pr 

^ [1 - 7(p, + Pr) - 7(pe2 + Pr') + (p/ + Pr') 
+PcPr(l2-7(pe + Pr))] 

3A 

^i-4[39 + 123pc - 97pr - 9pc' - 37pr' + 3(pe' + Pr') 



-pePr(l6 + 2l(pe + Pr))] 



3A 

+ -4[47 + 95pe - 61pr - 13pe' - Pr' + '^{p,'' + Pr' 

Ami 



+PcPr(32-2l(pe + Pr))] 
+{3pe'(l - Pr') 

3Ai 



2ml 



(2 + 26p, - 4p, + llpe' + 2pr' - 36pePr + HpcPr' + 9peVr') 
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xln 



6A2 
2 ' 

b 

(1 + - + VA)2 



+ ^(-1 - 8pc + 2p^ - 2pc^ - Pr"^ + 4pcPr + 2PcPt^ " SpcVr^)} 



4pc 



3A 

+{6p,2(l - Pc') + ^Pr'(46 - 28pe - 18pe') 



ISA 

f^p,2(l + 2pe-3pc')}ln 



(1 + p, -Pc + V^)^ 
4pr 



(C.l) 



with 



A = l-2(pe + p^) + (Pc-p^)' 



(C.2) 
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FIGURES 

FIG. 1. Relevant diagrams to obtain amplitudes T'^", T'^, T. (a), (b) correspond to W exchange 
terms, (c) - (f) the interference terms and (g), (h) charged Higgs exchange terms. 
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